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1  Abstract 


It  is  shown  that  the  direct  correlation  function  of  a  mixture  of  hard  ions 
in  the  mean  spherical  approximation  of  Hiroike  [lj,(MSA)  can  be  expressed  in 
terms  of  overlap  functions  of  charged  spherical  spherical  shells.  In  particular, 
if  the  system  has  a  mixture  of  pairs  of  ions  of  equal  size  and  opposite  charge, 
then  the  MSA  direct  correlation  function  is  given  by  the  electrostatic  energy 
of  a  pair  of  charged  shells,  of  radius  equal  to  the  radius  of  the  hard  ion  plus 
As  in  the  theory  of  Rosenfeld  [2]  this  direct  correlation  function  can  be 
derived  from  a  free  energy  functional,  and  a  simple  extension  to  non-uniform 
systems  is  given. 

2  Introduction 

Liquid  state  ‘heories  like  the  MSA  and  the  HNC  can  be  derived  as  varia¬ 
tional  problems  of  the  free  energy  functional,  which  is  written  in  terms  of 
the  Ornstein-Zernike  direct  correlation  function  .  This  view  was  recently 
introduced  by  one  of  us  [2,  3,  4],  for  systems  of  hard  objects  in  general,  and 
hard  spheres  in  particular  .  In  that  work  the  relation  between  the  grand 
potential,  the  direct  correlation  function  and  the  scaled  particle  theory  was 
clearly  shown.  We  extend  now  this  analysis  to  the  case  of  charged  hard 
spheres  in  general,  and  to  the  primitive  model  of  electrolytes  in  particular. 
The  asymptotic  limit  of  strong  Coulomb  interactions  between  the  charged 
particles  ,  that  is  the  limit  in  which  either  the  charge  goes  to  infinity  or 
the  temperature  goes  to  zero,  is  the  starting  point  of  our  present  discussion 
(5,  6]:  In  this  limit  the  free  energy  and  the  internal  diverge  to  the  same  order 
in  the  the  coupling  parameter  while  the  entropy  diverges  at  a  slower  rate 
and  therefore  the  free  energy  and  the  energy  coincide.  The  mean  spheri¬ 
cal  approximation  (MSA)  and  the  hypernetted  chain  approximation  (HNC) 
coincide.  This  is  a  very  gratifying  feature,  because  the  HNC,  which  from 
the  diagram  expansion  (and  numerous  test  cases)  point  of  view  is  the  more 
accurate  theory,  is  in  general  difficult  to  solve  numerically,  while  the  MSA 
is  analytical  in  most  cates,  and  in  the  asymptotic  limit,  of  a  rather  surpris¬ 
ingly  simple  form.  In  the  asymptotic  limit  the  excess  electrostatic  energy  is 
identical  to  the  exact  Onsager  lower  bound,  which  is  achieved  by  immersing 
the  entire  hard  core  system  in  an  infinite  neutral  and  perfectly  conducting 
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(liquid  metal)  fluid.  The  Onsager  process  of  introducing  the  infinite  con¬ 
ductor,  naturally  decouples  all  the  different  components  in  the  system  which 
may  differ  in  size,  shape,  charge  distribution  an  relative  orientation  in  space. 
As  a  result,  the  variational  free  energy  functional  in  the  high  coupling  limit 
diagonalizes,  and  the  mathematical  solution  of  the  asymptotic  problem  is 
given  in  terms  of  geometrical  properties  of  the  individual  particles  in  the 
system.  This  novel  approach  is  based  on  the  finding  that  both  the  weak  and 
strong  coupling  solutions  of  the  HNC  equation  are  bounded  from  below  by 
these  limits.  This  is  an  exact  new  result.  The  Onsager  approach  to  charged 
hard  particle  systems,  as  outlined  above,  has  been  developed  for  plasmas  of 
various  kinds  (charges  of  one  sign,  with  no  hard  core  excluded  volume  con¬ 
straints,  and  in  an  uniform  neutralizing  background  of  opposite  sign),  and 
for  uncharged  hard  particle  systems  both  in  the  uniform  an  non  uniform 
cases.  These  two  cases  provide  the  basis  for  the  treatment  of  the  more  gen¬ 
eral  cae  of  the  charged  hard  particle  system.  As  it  has  been  shown  in  the  case 
of  the  non-uniform  hard  particles,  this  new  approach  provides  an  excellent 
starting  point  for  the  quantitative  discussion  of  inhomogeneous  uncharged 
hard  sphere  systems.  We  show  that  the  MSA  free  energy  is  the  appropiate 
functional  from  which  the  charge  part  of  the  direct  correlation  function  is 
obtained  by  functional  differentiation.  This  also  suggests  a  simple  extension 
of  the  pair  MSA  to  inhomogeneous  charged  systems. 


3  Direct  Correlation  Functions  in  Terms  of 
Geometry  and  Electroststics 

Consider  an  arbitrary  mixture  of  charged  hard  spheres  of  radius  Ri  = 
er,/2,  charge  z,  and  number  density  p,  =  IVj/V.  The  temperature  is  T  and 
Boltzmann’s  constant  is  Jcb ■  We  use  the  notation  /3  =  We  can  imag¬ 
ine  that  the  charges  in  our  system  can  be  turned  on  gradually:  Then,  the 
Helmholtz  free  energy  F  is  the  sum  of  two  contributions,  a  term  due  to  the 
excuded  volume  effects  of  the  hard  spheres  Fhl  and  a  term  due  to  the  charges 
F ,ch,  which  usually  will  include  cross  terms  of  charge  and  excluded  volume 
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(cavity  effects). 


F  = 


F^»  _|_  piMSA 


(1) 


The  hard  sphere  term  was  discussed  in  detail  in  the  original  work  of  Rosenfeld 
[2,  3,  4],  and  can  be  written,  for  the  general  inhomogeneous  system 


Fmsa  =  J  dr i  $ch(r) 


(2) 


where  $ch(r)  is  the  local  Helmholtz  free  energy  per  unit  volume.  In  this  note 
we  use  the  mean  spherical  approximation  (MSA)  to  compute  this  magnitude: 
In  the  homogeneous  bulk  phase  [1,  8] 

$d>  ^  $MSA  _  gcharge  _|_  p3  j 3^  (3) 


where  the  local  internal  energy  per  unit  volume  is 

0Eehar*e  =  io'EpiZiNi 

i 

with 

_  e2  _  fle2 
tokfjT  eo 

In  this  equation  the  parameter  iV,  represents  on  half  of  the  potential  of  a 
spherical  shell  of  radius  6,  =  R,  +  1/2I\  This  energy  is  the  sum  of  the  energy 
of  a  collection  of  uncorrelated  spherical  capacitors  of  radius  We  could 
use  the  Onsager  picture  of  spherical  shells  in  liquid  metal  to  represent  this 
system. 


N-  =  ^ 

‘  i  +  r<7, 

(6) 

where  the  effective  charge  is 

A  =  2.  +  Wi 

(7) 

and  t]  is  defined  by 

r  l+rv, 

(2/t)(1  -  6)  +  Z)  T& 

(8) 

with 

6  =  (»/6)X>ft3 

(9) 

(4) 

(5) 
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The  parameter  rj  is  generally  small  for  ionic  solutions,  and  is  zero  for  the 
restricted  case  of  ions  of  only  one  diameter,  and  also  when  neutral  pairs  of 
ions  are  of  the  same  size.  The  paramenter  T  is  obtained  from  the  optimization 
condition 

d<bch 

^r  =  °  <10> 

which  turns  out  to  be  exactly  equivalent  to  the  equation  obtained  directly 
from  the  boundary  conditions  of  the  MSA  [7,  8].  The  energy  in  the  MSA 
can  be  considered  as  the  sum  of  the  energies  of  spherical  capacitors  of  ca¬ 
pacitance  C,  =  bx [9].  The  extension  of  4  to  non  uniform  systems  is  quite 
straightforward.  Consider  the  total  internal  energy  U 


UMSA  =  J  drlEMSA(rl) 

where  the  local  internal  energy  is 


(11) 


(12) 


We  rewrite  this  equation  using  the  charge  neutrality  sum  rule  [10,  11,  12]  for 
the  inhomogeneous  system 


-  =  J  pI(rt)zjh,i(ruT2) 


(13) 


Here  the  function  prevents  the  overlap  of  ions.  In  the  Onsagerian  picture 
of  the  MSA,  the  plausible  form  of  this  function  corresponds  to  completely 
uncorrelated  distribution  function  for  non  overlapping  pairs.  Therefore  this 
function  is  -1  when  the  distance  r  <  <T,y ,  and  zero  otherwise.  Substituting  13 
into  11  yields 


I 1U 


MSA 


=  -7o  I  dr\  I  dr2j2z*Njp,(ri)pj(r2)htJ(rur2)  (14) 

J  J  >,j 


Consider  now  the  interaction  of  two  charged  shells  of  radius  bj  ,  with 
centers  at  positions  ri,r2.  From  simple  electrostatics  we  know  that 

bj ;  r,2)  =  (1/2)  J  dr3?,(ri  -  r3)^(r3  -  r2)  r,2  =  rx  -  r2  (15) 
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where  the  electrostatic  potential  4>  is  defined  by 

I  dr>r  -i  (i«) 

j  I  r32  —  ^*4  | 

The  charge  density  of  the  shells  is 

*<r>  =  J^(l  r  1  -b')  (17> 

The  total  energy  is  given  by  14 

0UMSA  =  -(l/2)7o / J dr2'Epi(Tl)pi(T2)9ij(bi,bj-,rl2)  (18) 


where  we  have  made  use  of  the  extension  of  the  electroneutrality  relation  13 
to  the  case  of  smeared  charge  distributions.  Consider  again  the  interaction 
of  the  spherical  shells  15:  When  they  overlap  we  get 


*«(».,  h-,  r)  =  -  b,f  -  2r(b,  +  b,)  + 

0,0 jr 

(19) 

If  we  rearrange  this  expression,  we  obtain 

PIjmsa  =  _(i/2)70  /  dr3£Q,(T3)i ,(r3) 

•J 

(20) 

where 

Qi( r3)  =  J  drlpi(ri)q,(r1  -  r3) 

(21) 

and 

$j(r3)  =  J  dr2pj{r2)<i>j(r3  -  r2) 

(22) 

For  ryz^jrxiy©da^)hig^eQn4guia4itMisva*(&kspatially  uniform  systems  20  is  just 
equal  to  the  starting  equation  4.  In  fact,  simple  integration  of  2fand  22  with 
17  yields 

Q,(  r3)  =  (23) 

and 

*i(r3)  =  -2  Nj  (24) 
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We  recall  now  the  relation  between  the  Helm^holtz  free  energy  and  the 
direct  correlation  function 


rMSA 


(rx,r2)  -  - 


62(3F  MSA 

^P.(r  i)Spj(r2) 


(25) 


Let  us  now  restrict  to  the  case  in  which  r)  —  0. Functional  differentiation  will 
yield 

r>)  =  (l/2h»*ij  (26) 

which  is  exactly  identical  to  the  direct  correlation  function  first  derived  by 
Hiroike  [1] 


Equation  26  is  the  natural  extension  of  the  MSA  direct  correlation 
function  to  inhomogeneous  sy terns.  To  a  first  approximation  we  may  take 
the  bulk  value  of  1',  although  refinements  are  possible.  The  pair  correlation 
functions  which  are  computed  from  the  Ornstein  Zernike  equation,  can  be 
shown  to  satisfy  the  local  electroneutrality  condition  13.  A  more  detailed 
discussion  will  be  published  in  the  near  future. 
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